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Abstract. We probe the U(N) chiral Gross-Neveu model with a source- 
term J\E n I'. We find an expression for the renormalization scheme and scale 
invariant source J, as a function of the generated mass gap. The expansion of 
this function is organized in such a way that all scheme and scale dependence 
is reduced to one single parameter d. We obtain a non-perturbative mass 
gap as the solution of J = 0. A physical choice for d gives good results 
for N > 2. The self-consistent minimal sensitivity condition gives a slight 
improvement. 
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1 Introduction 



In a previous paper, we developed a method for dynamical mass generation 
in asymptotically free quantum field theories. It was applied to the ordi- 
nary Gross-Neveu model, Jjj], and a mass gap was found, Q, which agreed 
very well with the exact result, ||. In this paper we will apply the same 
method to the non-abelian Thirring model (NATM) or chiral Gross-Neveu 
model (CGNM), [ffl. It is another one of those rare quantum field theories 
where exact results, like the mass gap, can be obtained. In Q the mass 
gap is calculated exactly in terms of A which is the non-perturbative mass 
parameter which sets the scale for the running coupling in a certain scheme. 
Comparing our results with j|] will provide another check on the accuracy 
of our method. 

The idea behind the method is very simple. A source term, J^^f, is 
added to the NATM-Lagrangian and then we calculate the mass gap using 
ordinary perturbation theory to obtain the perturbative expansion for m(J). 
As a consequence of asymptotic freedom, this expansion is only valid for 
large values of J. If we let J approach zero, the coupling constant grows too 
large, and the perturbative expansion for m(J) becomes invalid. Therefore, 
we cannot take the limit J — > 0. If instead we consider the perturbative 
expansion for the inverted relation J(m), perturbation theory remains valid 
in the limit J —* 0, provided that a solution m exists for J{m) = 0, which is 
not too small. As in ordinary perturbation theory, the result for the mass 
gap m is renormalization scheme (RS) and scale dependent. To eliminate the 
mass renormalization dependence we use the scheme and scale independent 
quantity J, instead of J. Exchanging g 2 (fi) for l/(/?o hi a 7 ) as the expansion 
parameter, reduces the remaining dependence to one single number d, which 
can be fixed by some external physical condition, or in a more self-consistent 
approach, by the principle of minimal sensitivity. 

The paper is organized as follows. In section [2] we will present the results 
necessary for application of the source inversion for the NATM. In order to 
avoid unnecessary repetition, we shall refer to the paper Q for the derivation 
of the general formula. The outcome of our calculations will be discussed 
in section ||. As a bonus we will show that reparametrization of the d- 
dependence will enable us to solve the mass gap equation exactly. Details 
of the exact evaluation of the finite parts of the two loop Feynman integrals 
which occur in the sunset topology are given in an appendix. 
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2 The non-abelian Thirring model 



The U (N) invariant NATM describes the interaction of iV single flavor Dirac 
fermions Vl/ a , a = 1, . . . , N in two dimensions with the (massless) Lagrangian 



(2.1) 



where T\ i = 1, . . . , N 2 - 1, are the generators of SU(N) with the normal- 
ization Tr(r i T J ') = \ ■ (Note that our coupling constant g 2 is two times 
the coupling constant g 2 of ||.) This model is also known as the CGNM be- 
cause a Fierz-transformation of the interaction term leads to the equivalent 
Lagrangian 
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The NATM is asymptotically free, || and possesses, apart from the U(N) 
invariance, a chiral £7(1) symmetry. In ordinary perturbation theory this 
symmetry remains unbroken and no mass gap is generated. 

We begin by perturbing ( |2,1| ) with a fyfy composite operator to produce 
the new Lagrangian 



(2.3) 



The detailed three loop renormalization of this model using dimensional 
regularization has been given in || |?J which built on the one and two loop 
calculations of || !> H> tl- The results for the (3 and 7- functions in the 
MS'-scheme are (g 2 = g 2 (p)): 
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To apply the source inversion at two loop order we also require the two loop 
perturbative result for the mass gap. To determine this we have computed 
the fermion two point-function at two loops and extracted the finite part 
exactly after performing the renormalization. The values for the integrals 
we obtained have been checked against the numerical results of [10] for the 
mass gap of the ordinary Gross Neveu model. We find 
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where C(^) is the Riemann zeta function. From this one easily arrives at the 
expansion for the inverted relation 
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We define Xq and Yq as the coefficients which do not multiply a logarithm 

J(m) = m(l+ 5 2 (m)X + (7 4 (m)yo + ...) • (2-8) 

The expansion (2.7) is highly scheme and scale dependent. This dependence 
is reduced drastically if we replace J{n) with J which is the scheme and 
scale independent quantity associated with J, and then expand in powers of 
1/ ^/5oln^-^ rather than in g 2 {^)- Starting with the expansion for J(m) in 
a general scheme, we found [0] 
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All the scheme and scale dependence now resides in d = (3q In ( -j?r-^z) and 
we can recover the original NATM, by putting the naked source Jq equal to 



zero 



Jo(m) ~ J(m) = . 
We find a non-perturbative mass gap which is a solution of 
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The total series is of course (i-independent but one can only calculate it up 
to a certain order in perturbation theory which will give us a mass gap that 
depends on d. One can check that the (i-dependence of the order n truncated 

series is 0[ 3 . We will consider two possible ways of fixing d. 

V O In ~y hd / 

MS 

The first one reduces to a choice for A, that corresponds to a physical scheme. 
The second one fixes d by the principle of minimal sensitivity. In Q we used 

the value of the expansion parameter 1/ ( (5q In + d J as a source of error 

V MS J 

estimation. This works if the coefficients are of order one. Assuming that 
the series is asymptotic, a rather large value of the expansion parameter can 
still give reasonable results, as long as the complete terms in the series are 
small. In the next section we will show that this is indeed the case. For the 
2-loop results it is better to estimate the error from the second order term, 
than from the expansion parameter. 



3 Numerical results 

The exact result for the mass gap was obtained in 



e 2JV 



m 



iT Apy ' 



(3.13) 



where Apy is defined as the scale parameter for the running coupling, with 
a condition on the normalized four point function, calculated with a Pauli- 
Villars regularization. To obtain m/A-rj-j we need to determine the relation- 
ship between the renormalized coupling g of the dimensional regularization 
MS'-scheme and the coupling gpy used in the Pauli-Villars scheme. This 
can be achieved by comparing the normalized fermion four-point function 
to one loop order in both schemes. We find 
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(See, for example, plfl .) So we finally arrive at 
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where 7 is the Euler-Mascheroni constant. 
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Table 1: Physical scheme results. 



3.1 Physical scheme 

As in ||], we will define a physical RS on the normalized 4-point function of 
Cj, to obtain a physical value df for d. Demanding that g 2 coincides with 
the 4-point function at zero external momentum (see equation (5.5) of |Q]), 
one arrives at 

° 2 = 5 K i+ ^ + ---) (3 - i7) 

Taking fi 2 = m? leads to df = ^--[] We now find the one and two loop mass 
gaps mfi and m/2 as the solutions of the one and two loop truncation of 
( |2.12| ) with d = df. The deviations from the exact result for mji and my 2 
for the N — > 00 limit and the first order result in 1 /N have been displayed in 
table in terms of a percentage. We also provided the value of the two loop 
expansion parameter 1 / (2/3q In + df) = y and the second order term II 

MS 

of the series in ( |2.12j ). For N = 2,3 we find no one-loop mass gap. All the 
other results lie somewhere between the iV — > 00 and \/N approximations. 
We also observe convergence. In other words the comparison with the exact 
result improves for the two loop truncation. From the y/ir and the II columns 
we learn that II clearly gives a better indication on the size of the error. 



*As in |2|, one can show that every physical value for d (~ N) gives the correct N — > 00 
limit. 
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3.2 Minimal sensitivity 



The equation for m(d) ( [2.12 ) can only be solved numerically. If we consider 
instead the expansion parameter y as the free parameter one can solve it 
analytically to find m(y). Indeed, we can rewrite ( 2,12j ) as 
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The one- and two-loop truncation of (3.18) is now solved easily. At one loop 
it is a linear equation in k and one finds k = ^~(^ + ^4o)- After substituting 
this into (3.19) we find the one loop mass gap to be 
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The 2-loop truncation gives a quadratic equation in k, with two roots 
ki(y), k2(y)- Hence, the two solutions for the mass gap are 
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The behavior of mi(y) is more or less the same for all values of N. One ob- 
serves a sharp maximum, followed by an asymptotic descent to zero. There 
is no region of minimal sensitivity. For N > 2 the situation changes at 
two loops. One of the two solutions rri2i{y) has, in addition to the sharp 
maximum, a rather flat minimum. This is the point of minimal sensitivity. 
In figures [l] and |2| we plot the one- and two-loop solutions for the generic 
N = 5 case. 

The other 2-loop solution is not physical since it varies enormously in 
the region of interest, defined as the region with acceptable estimated error, 
and no minimal sensitivity is found. For N = 2 the two-loop solution has no 
minimum with instead only a rather sharp maximum at 68% deviation. No 
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Table 2: Minimal sensitivity results. 



true minimal sensitivity point can be identified. The results for N > 2 are 
displayed in table [T^. They are slightly better than the two-loop physical 
scheme. Again we find II to provide a better indication on the error then 
y/ir. We finally remark that also the minimal sensitivity condition can be 
solved exactly, to give an analytic form of the 2- loop mass gap. We will not 
present it here, however, since it is a large expression and does not give any 
new insights. 

4 Conclusions 

We have successfully applied the source inversion method to the chiral Gross- 
Neveu model. This required a two-loop calculation of the mass gap in the 
massive NATM which we carried out exactly. Comparison with the exact 
result for the non-perturbative mass gap gives a satisfying match. For the 
physical scheme, there is convergence of the 2-loop result versus the 1-loop 
result. The 2-loop results are good for N > 2 with a 10%-deviation for 
N=3 and < 5% for A > 3. The minimal sensitivity condition gives a slight 
improvement. As in the case of the ordinary Gross-Neveu model, the N = 2 
result is poor. The two-loop physical scheme gives a 46% deviation. The 
success/failure of the method for A > 2/N = 2 is fairly consistent with the 
error estimation one obtains from the second order term in the mass gap 
equation. Finally, it would be worthwhile to apply the technique discussed 
here to other models where exact mass gaps are also available. This would 
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have the long term aim of applying the procedure to theories where the only 
information on the dynamical generated mass comes from say Schwinger 
Dyson or lattice methods in order to ascertain how competitive the results 
would be. 
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A Computation of two loop integrals. 

In this appendix we discuss the evaluation of the basic Feynman integrals 
which underly the exact value of our mass gap at two loops. At one loop 
there is only one basic integral which is defined by 

1 = i I \U2 1 21 ( A ' X ) 

Jk [k z - m z \ 

in Minkowski space where Ju = / k / (2tt) 2 and it has the exact value in 
w-dimensions 

_ r(i- r /2),„<^ 

(4tt) w / 2 V ; 

Therefore, if uo = 2 — e then / has a simple pole in e which is the foundation 
of the one loop renormalization. At two loops all contributions to the 2- 
point function can be reduced to several basic Feynman integrals. These are 
I 2 , A(p 2 ) and A^(p 2 ) where 

„2\ -2 I" I" - 
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and 

Ik Ji [(k — p) 2 — m 2 ] [I 2 — m 2 ] [(k — I) 2 — m 2 ] 

and these latter functions only occur in the sunset topology. Other integrals 
with an obvious definition such as A^(p) and A fiucr (p) also arise but the 



relevant 2-point function contributions can be related to (|A.3|) and (A. 4) 



[12|. For instance, 

A^Cp) = 2 -p,A(p 2 ) , p"A^(p) = 2p^A, v {p 2 ) - 2 -p 2 {p 2 -m 2 )A(p 2 ) . 

(A.5) 
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It is elementary to observe that A(p 2 ) is finite in two dimensions. Hence, for 
the mass gap we only need to evaluate it in two dimensions when p 2 = m 2 . 
To do this we follow the Feynman parameter approach of |l4|] which gives 
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in kj-dimensions after carrying out the momentum integrations. Restrict- 
ing to two dimensions the y-integration can be performed from an integral 
representation of the hyper geometric function, 2-^1(^1 b; c; z), giving 
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Next we set p 2 = m 2 in the two dimensional integral to obtain 
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The final integral can now be calculated exactly, ||15|| , to produce 
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The remaining integral which arises in the sunset topology occurs with 
two Lorentz contractions. First, in ^-dimensions without setting the on-shell 
condition it is straightforward to show that, [12], 
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Although the contraction of (A. 4) with p^p v is also divergent it cannot be 
written in a similar closed form. However, its divergent part is known to be 
P 2 I 2 /uj, |2|. Therefore, 



P 2 K(P 2 ) 



(A.12) 
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will be finite in two dimensions and can be evaluated exactly when the on- 
shell condition is set similar to the derivation of ( A.lOj ). With the same 



Feynman parametrization as (A.6) we have 
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Hence, using the properties of the hypergeometric function again we find 

F A (P 2 ) , 2 2 = ft f'dx { -^\lnx (A.14) 
uj=2 , p 2 =m 2 8ir z Jo (1 + X) 



leading to, [15| , 



Fa(p 2 ) 9 2 2 = £[1 - C(2)] • (A.15) 

Hence, all integrals in the full 2-point function can be written in terms of 
I 2 , flOcD and ( |Q5D . 



We have checked that the values obtained here for the finite parts of 
A(p 2 ) and A^ u (p 2 ) agree with the numerical values given in |T^|. For in- 



stance, repeating the calculation which leads to equation (4.5) of |1C] we 
find that [0.737775 - vr 2 /96] corresponds to [C(2)/2 - 3/16]. 
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